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Abstract. We present the dynamical structure of the phase space of 
p ' the planar planetary 2/1 mean- motion resonance (MMR). Inside the 

r-["1 | resonant domain, there exist two families of periodic orbits, one asso- 

ciated to the librational motion of the critical angle (c-family) and the 
other related to the circulatory motion of the angle between the peri- 
centres (Aro-family). The well-known apsidal corotation resonances 
(ACR) appear at the intersections of these families. A complex web of 
£^ ■ secondary resonances exists also for low eccentricities, whose strengths 

and positions are dependent on the individual masses and spatial scale 
of the system. 

Depending on initial conditions, a resonant system is found in one 

■ of the two topologically different states, referred to as internal and 

external resonances. The internal resonance is characterized by sym- 

' metric ACR and its resonant angle is 2 A2 — Ai — w\ , where Aj and 

mi stand for the planetary mean longitudes and longitudes of pericen- 
tre, respectively. In contrast, the external resonance is characterized 
by asymmetric ACR and the resonant angle is 2 A2 — Ai — -072. We 
show that systems with more massive outer planets always envolve in- 
side internal resonances. The limit case is the well-known asteroidal 
resonances with Jupiter. At variance, systems with more massive in- 
ner planets may evolve in either internal or external resonances; the 
K> \ internal resonances are typical for low-to-moderate eccentricity config- 

■ urations, whereas the external ones for high eccentricity configurations 

of the systems. In the limit case, analogous to Kuiper belt objects 
in resonances with Neptune, the systems are always in the external 
resonances characterized by asymmetric equilibria. 



1. Introduction 

Asteroids and exoplanets are amongst the most striking sources of problems in 
Celestial Mechanics where chaos and order play an important role. Asteroids 
have always puzzled astronomers by their peculiar distribution: in some regions, 
they avoid resonances, but in others, they prefer to be exactly where the reso- 
nances are found. The famous problem of the Kirkwood gaps, solved in the last 
quarter of the past century, is the better example. Asteroids have benefited from 
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more than one century accurate observations and the number of asteroids whose 
orbits are known with high precision is presently of order of many thousands. Ac- 
curate observations allowed us to map the phase space where the asteroids evolve 
and obtain the boundaries of regular and chaotic motions. In addition, the in- 
creasing computer capacity allowed us to understand why they are distributed 
in such a way. 

Exoplanets on the other hand are being observed at the limit of our technical 
capabilities. We are far from having accurate knowledge of their orbits. We know 
presently more than 700 planets. Most of them are in systems composed of only 
one planet. But the number of multi-planet systems is increasing and may reach 
its first hundred in a short time. Besides, what is important, they show a striking 
diversity. They range from the small systems of super-Earths and mini-Neptunes, 
as those discovered by the space missions CoRoT and Kepler in circular orbits 
very close to their central stars, to systems of giant planets in very elongated 
elliptic orbits reaching enormous distances from their host stars, as far as the 
outskirts of our Solar System. Eccentricities larger than 0.5 are not infrequent: 
in multi-planet systems they reach up to 0.75 and in single-planet systems they 
reach up to 0.97. 

Many multi-planet extrasolar systems exhibit resonant behavior. Systems 
of resonant planets is a novelty for the Celestial Mechanician. We were used 
to study the resonant problems of small objects, whose mass can be assumed 
as zero (restricted three-body problem), or a few planetary satellites in which 
restricted models cannot be used but which are dominated by the oblateness of 
the central planet. So, exoplanets are a prime source for cases of the three-body 
problem in which we cannot neglect the mass of the bodies and in which the 
mutual point-like gravitational interaction is by far the dominant force. The 
most frequent resonance found is the 2/1 mean-motion resonance (e.g. GJ876 
c-b, HD 40307 c-d, HD 73526 b-c, HD 82943 c-b and HD 128311 b-c). 

Although other resonance sites may also be inhabited (e.g. HD 60532 b-c 
in the 3/1 or HD 45364 b-c in the 3/2 resonance), the 2/1 is presently the most 
populated. Notwithstanding the limited accuracy of the observations, resonant 
configurations have at least two strong arguments in their favor: The long-term 
stable motion of close planets in high eccentricity configurations is possible only 
if the planets are locked in (and protected by) MMRs. Secondly, resonance 
trapping appears to be a natural outcome of planetary migration processes due 
to planet-disk interactions, which are believed to take place in the latest stage 
of the planet system formation (e.g. Kley 2000; Snellgrove et al. 2001; Kley et 
al. 2005). 

The classical studies of the general three-body problem with their aver- 
aged Hamiltonians have been quickly extended to the case of exoplanets (Had- 
jidemetriou 2002; Lee & Peale 2002; Beauge k Michtchenko 2003; Ferraz-Mello 
et al. 2003). The periodic solutions, dubbed as corotation resonances or apsi- 
dal corotation resonances (ACR), were obtained to cover all possible families in 
low-order resonances (Lee 2004; Beauge et al. 2006; Michtchenko et al. 2006b, 
Giuppone et al. 2010). ACR are stationary states of the resonant Hamiltonian 
averaged over the synodic period and correspond to equilibrium solutions of the 
averaged equations of motion. 
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Notwithstanding the attention devoted to ACR and periodic orbits, not 
much was known on the topology of the phase space outside their vicinity. A 
detailed analysis was presented first by Michtchenko & Ferraz-Mello (2001) for 
the 5/2 MMR, by Callegari et al. (2004) for the 2/1 MMR and by Callegari et al. 
(2006) for the 3/2 MMR. However, these works were limited to the study of a few 
specific systems and are only valid for very small eccentricities. Even so, results 
showed the complex structure of the phase space, populated by several different 
families of periodic orbits, modes of oscillation and possible regimes of motion. 
ACR appear to be only one of several distinct types of stable configurations. 

It was necessary to extend our knowledge of the resonant dynamics of the 
MMR beyond the ACR. This has been done for the 2/1 resonance in two papers 
published by Michtchenko et al. (2008 a,b), for arbitrary mass ratios and with 
no restrictions in the orbital eccentricities, providing a complete map of the 
phase space around the ACRs. These results, summarized in this paper, help us 
understand the results of numerical simulations and pinpoint where future (or 
even currently known) extrasolar systems may be found. Even if all planetary 
systems in this commensurability seem to be restricted to the very close vicinity 
of ACR, the maps allow to see the possible evolutionary routes inside the 2/1 
MMR from initially non-resonant configurations. 



2. The model of the resonant three-body problem 

When the ratio of the orbital periods of two planets is close to a ratio of two 
simple integers, we say that the planets are close to a mean-motion resonance 
(MMR). A mean- motion resonance is often written in a generic form as 



m/n 2 = {p + q)/q, 



where nj are the planetary mean motions and p and q are integers, the latter one 
representing the order of the resonance. The critical (or resonant) angles 0j are 
defined as 

01 = {l+p/q)X 2 -{p/q)Xi-w 1 , ,^ 

02 = (1 +p/q)M - (p/?)Ai - w 2 , K > 

where \ and Wi are mean longitudes and longitudes of pericentre of the planets; 
index 1 denotes the inner body, while index 2 is reserved for the outer body. 
The behavior of the critical angles define the location of the system with respect 
to the resonance: when one of these angles is librating, the system is said to 
be inside the resonance. It is worth noting that the secular angle defined as a 
difference in longitudes of pericentre of the planets, is Azu = o\— cr 2 . 

The dynamics of two resonant planets, with masses mi and m 2 , orbiting 
a star of the mass mo is defined by the averaged Hamiltonian Tires and two 
integrals of motion, AM and /C, where the first is the total angular momentum 
and the second is the so-called spacing parameter. These analytical functions 
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are given by the expressions, up to second order in masses: 

%rcs = ~^2i=l 2 cii 1 ~ 2tt So ^( a it e ii °l) Q) dQ , 



AM = mini a? Jl - e\ + m 2 n 2 a 2 \J 1 - e\ ( 2 ) 

K. = (p + q) mini a\ + p m 2 n 2 a2, 

where 1Z is the disturbing function. All orbital elements, including the semi- 
major axes cjj and eccentricities ej, are canonical astrocentric variables; their 
relation to the usual osculating heliocentric orbital elements can be found in 
(Ferraz-Mello et al. 2006). The averaging of the disturbing function 1Z is done 
with respect to the synodic angle Q = A 2 — A] . In the vicinity of the MMR, 
Q is much faster than the resonant and secular angles, and does not influence 
significantly the long-term evolution of the system. Thus, all periodic terms de- 
pendent on Q can be eliminated (i.e. averaged out) of the Hamiltonian function, 
and only secular and resonant terms need to be retained. 
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Figure 1. Numerical simulation of two planets displaying large- 
amplitude librations around a (0, 0)-ACR, with masses and initial con- 
ditions drawn from Fit B of Ferraz-Mello et al. (2005) for the HD 82943 
system. Plot shows the resonant evolution of the semimajor axes (left 
panel) and eccentricities (right panel), averaged over short-term oscil- 
lations. It is worth comparing the resonant oscillation timescale (left 
panel) to the secular one (right panel). 

The invariance of AM and K, has important consequences for the orbital 
evolution of the system. It indicates that, after the averaging process, from 
the set of the variables (ai, o 2 , ei, e 2 ), only two are independent, and the pla- 
nar resonant problem has two degrees of freedom. Each degree of freedom is 
characterized by its proper mode of motion. Note that K, depends only upon 
masses and semimajor axes; therefore, the invariance of the spacing parameter 
K, defines a coupling of the semimajor axes: they oscillate with opposite phases 
and with amplitudes that are inversely proportional to the planetary mass. We 
say that the semimajor axes have no secular variation and their evolution occurs 
in timescales associated with the resonant mode of motion. This effect can be 
seen in Figure Q] (left panel) . 

On the other hand, the integral of total angular momentum AM causes a 
coupling of the planetary eccentricities, in such a way that, when one eccentric- 
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ity grows, the other decreases, as seen in Figure [T] (right panel). The oscillation 
generally occurs in timescales much longer when compared to those of the reso- 
nant mode. This mode is defined as secular mode of motion of the system and 
is associated with the secular angle Aw. AM is a function of both semimajor 
axes and eccentricities, that implies that both modes of motion, resonant and 
secular, are present in eccentricity variations. In the example shown on the right 
panel in Figure [IJ the amplitudes of the resonant component in the eccentricity 
variations are much smaller (and frequencies are much higher) than those of the 
secular mode. However, as will be shown below, for some initial conditions, both 
amplitudes and frequencies can be comparable. 

3. Stationary orbits of the averaged 2/1 resonant problem 

The Hamiltonian function given in Eq.([2]) is very complicated, even in the planar 
case, and generally possesses several extrema, which define stationary solutions of 
the averaged problem. These solutions are often referred to as Apsidal Corotation 
Resonances (ACR). To obtain these special solutions over a large domain of 
the parameters of the problem, we employ the geometrical method presented in 
Michtchenko et al. (2006b), where each equilibrium solution is identified looking 
at local maxima of the Hamiltonian function, for given values of the total angular 
momentum AM and the spacial parameter 1C. The location of the global maxima 
of the 2/1 resonance Hamiltonian calculated for all possible values of AM and 
a fixed /C, are shown on the plane (ex, e<i) in Figure [2j in the form of families 
parameterized by the mass ratio mijm\. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 00 02 °" 

e ! inner planet eccentricity 



Figure 2. The 2/1 resonance stable apsidal co-rotation solutions pa- 
rameterized by the mass ratio mj/m]. The cases with m^j 'mi > 1 are 
shown on the left panel, while with mj/mi < 1 on the right panel. The 
mass ratio 1.064 corresponding to the HD 82943 system is illustrated 
by a gray curve. 

Usually, ACR are classified in two types: symmetric and asymmetric solu- 
tions (e.g. Beauge et al. 2003, Lee 2004, Voyatzis and Hadjidemetriou 2005). 
The symmetric solutions are characterized by stationary values of both criti- 
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cal angles ([T]) at or 180°. The symmetric ACR-solutions evolve monotonically 
with the increasing eccentricities of the planets. For instance, for milm\ > 1.015 
(Ferraz-Mello et al. 2003), all ACR families are symmetric; some are plotted on 
the left graph in Figure [2J 

For mass ratios smaller than this critical value, the smooth evolution along 
a symmetric family is interrupted by a sudden increase of the outer planet's 
eccentricity (the case of the family with m2/m\ = 1 on the left graph). The 
corresponding symmetric solutions become unstable and the stationary systems 
evolve now along an asymmetric segment of the corresponding family. Asymmet- 
ric ACR are characterized by stationary values of both critical angles different 
from zero or 180°. The families shown on the right graph in Figure [2] exhibit 
this feature. 

A second critical value of the mass ratio occurs at ~ 0.36 and is illustrated 
by cyan and magenta curves in Figure [2] (right panel). Below this limit, asym- 
metric ACR-families split into two disconnected branches on the (ei, e2)— plane, 
marked by labels A and B. In these cases, there may exist two distinct stationary 
configurations described by different sets of orbital elements, but leading to the 
same value of total angular momentum. 

We analyze the evolution of the proper periods (inverse of the proper fre- 
quencies) along each ACR-family, employing a method of dynamic power spectra 
(Michtchenko et al. 2006a). The left graph in Figure [3] shows, in logarithmic 
scale, the period of the resonant proper mode of motion (T CT ) and the period of 
the secular mode (Ta w )> both obtained along the ACR-family parameterized by 
the mass ratio of the HD 82943 system. The ratio between both periods, plotted 
by a dashed line on the same panel, indicates that the characteristic times of the 
resonant component of motion are generally shorter (at least in order 1) than 
the secular ones. Both periods are comparable solely for very small eccentricities 
when, for some initial configurations of the planets, they can be commensurable 
originating a complex structure of secondary resonances inside the primary 2/1 
resonance. 




Figure 3. Left: Evolution of the period T a of the resonant mode of 
motion and the period T^vj of the secular mode along the ACR family 
parameterized by mv/mi = 1.064, with the individual planet masses 
of the HD 82943 system. The ratio T^ w jT G is shown by dashed line. 
Right: Evolution of the secular period 7a ro along several ACR families. 
The values of the mass ratio are indicated by numbers. 
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The right graph in Figure [3] shows the secular period, in logarithmic scale, as 
a function of the outer planet eccentricity obtained along several ACR-families 
with different mass ratios. For large values of 777,2/7711, there is a little change in 
the period with a maximum value at e<i roughly between 0.1 and 0.2. However, 
for decreasing values of 777,2 /mi, the maximum becomes more pronounced until 
it reaches a singularity near e2 — 0.1, for the mass ratio 771,2/mi = 0.015. In 
Figure [3] (right panel) we show discontinuities which appear along the ACR- 
family parameterized by 7772/7771 = 1. A comparison with the plot of Figure [2] 
shows that these discontinuities are associated to the bifurcation of the symmetric 
ACR leading to the origin of the asymmetric stable solutions. The existence of 
the singularities indicate that the difference between symmetric and asymmetric 
ACR is not restricted to the equilibrium values of the angles, but constitutes 
qualitatively distinct solutions, each belonging to separate regions of the phase 
space. 




Figure 4. Examples of stable symmetric (a) and asymmetric (b) and 
(c) stationary configurations of two planets in the 2/1 resonance, with 
the mass ratio 771,2/mi = 0.3. The continuous lines are apsidal lines, 
while the dashed lines are conjunction lines (in the case (a) two lines 
coincide) . The full circles show the positions of the planets on the orbits 
at the instant of the closest approach between them. 



4. Geometry of the stationary orbits 

Figure [J] shows three different configurations of the ACR-solutions of the 2/1 
MMR, obtained for 777,2/771,1 = 0.3. The case (a) illustrates a symmetric solution 
with anti-aligned pericentre lines. In this case, the closest approach between 
two planets occurs when the planets are in conjunction: the inner planet in the 
pericentre and the outer planet in the apocentre of their orbits (Aw = 180°). 
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Two asymmetric ACR-solutions are shown in Figure H](b) and (c). Both 
configurations are characterized by the same set of the constants AM and IC, but 
the solution (b) belongs to the branch A of the ACR-family, while the solution 
(c) corresponds to the branch B in Figure [2] (right panel) . The secular angle Aw 
defines the relative orientation of the pericentre lines of two orbits. The positions 
of the planets on their orbits are strongly tied by the 2/1 resonance relationship 
between the planetary mean anomalies which are related to the critical angles 
through Mi = o\ + (p + q) Q and Mi = &2 + pQ,where Q is the synodic angle 
and p = q = 1 (Beauge and Michtchenko 2003). Full circles in Figure [4] show 
positions of the planets, for which their mutual distance is minimal. 

In the configuration (b), when the orbits are not intersecting, the closest 
approaches occur when the planets are in conjunction, forming a co-linear con- 
figuration with the central star. It is interesting to observe that the planetary 
conjunctions avoid positions between the two pericentra. As a consequence, con- 
junctions always occur when two planets are far away from positions where the 
two orbits are very close one to another; this feature can be seen as a protection 
mechanism due to the 2/1 mean- motion resonance. 

In the case (c), two orbits are intersecting and the mutual distance is not 
more minimal when the planets (open circles) are in conjunction. The closest 
approaches occur when the central star and the two planets (full circles) are in a 
triangle configuration and the outer planet is localized close to the intersection of 
the two orbits. In this case, again, the protection mechanism prevents very close 
approaches between two planets (see also Lee et al. 2006, for similar analysis). 



5. Two regimes of resonant behavior: interior and exterior reso- 
nances 

In the previous sections, we have seen that the resonant behavior is actually 
composed of two proper modes of motion, the resonant and secular ones. It 
is generally characterized by the libration of one angle and the circulation of 
the other. We refer to the librating angle as the resonant angle and to the 
oscillation/circulation angle as the secular angle. The secular theories provide us 
with the secular angle: it is defined by the difference in longitudes of pericentre 
Aw = vo<i — w\ or, from our definitions ([1]), as Aw = a\—a2- The last expression 
shows that only one of two critical angles defined in Eq. [I] is independent and is 
a truly resonant angular variable of the problem (the other independent angular 
variable is Aw). The question which rises now is which one of two critical angles 
can be chosen resonant variable. 

When the outer body is more massive (mi/mi > 1.015), the situation is 
similar to that of an asteroid evolving in a resonance with Jupiter, when the 
librating angle is o~\. In this case, all ACR-solutions are symmetric. Thus, we 
choose the critical angle o~\ as a resonant angular variable of the problem and 
adopt the name interior resonance for this regime of motion. 

For the planets with the mass ratios 777.2/7711 < 1.015, both kinds of ACR- 
solutions are possible, symmetric and asymmetric ones. If 7712/7711 << 1, which 
is the case of a Kuiper belt object in resonance with Neptune, the asymmetric 
solutions are dominating. Therefore, in the case of asymmetric ACR, we choose 
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the critical angle a% as a resonant angular variable of the problem. We refer to 
the regime of motion characterized by asymmetric ACR as exterior resonance. 

The behavior of the secular angle Aw has also some special characteristics. 
In both interior and exterior resonances, it may be either a circulation or an 
oscillation. However, the oscillation in this case should be differentiated from 
the regime of motion defined as a libration in Celestial Mechanics. Indeed, the 
motions are akin to a family of concentric curves around one center displaced 
from the origin; the curves which are close to the center, do not enclose the 
origin - they correspond to oscillations, while the outer curves enclose the origin 
and correspond to circulations. The separation between them is not a dynamical 
separatrix, but just one curve passing through the origin. The whole set of 
curves forms a homeomorphic family of solutions and the distinction between 
oscillations and circulations in this case is merely kinematical. To stress this 
behavior we use the word oscillatiory /circulatory, when describing this regime 
of motion. 



6. Dynamics around the ACR solutions: Interior resonance 

Independent of the mass ratio, the dynamics associated to all symmetric ACR 
is qualitatively similar, showing interior resonances. To portray the phase space 
of the interior resonance, we introduce a symmetric representative plane {n\jni, 
e^), with the resonant angle a\ fixed at and the secular angle Aw at (positive 
values on the e2-axis) or 180° (negative values on the e2-axis). A typical example 
of the dynamical map of the phase space around one symmetric ACR is shown 
in Figure [5j 

The map contains one center, which represents the ACR solution; its posi- 
tion shown by a red circle. Around the center, we find a region of quasi-periodic 
motion, coded in gray scale on the dynamical maps. Light gray color indicates 
regular motion, while darker tones correspond to increasingly chaotic motions. 

The quasi-periodic behavior is the composition of two independent modes 
of motion: the resonant (with the frequency / CT ) and the secular (with the fre- 
quency /Aro) ones. The first one is associated to the resonant angle o~i, while 
the second one to the difference in longitudes of pericentre, Aw. Thus, any 
regular solution will be given by a combination of two periodic terms and their 
harmonics, each with a given amplitude and phase angle. Generally, both modes 
are well separated in the frequency space, with f a much higher than the secular 
frequency /a to (see Figure [3] left panel). However, at small eccentricities, both 
frequencies may have same order, and it is possible to find initial conditions cor- 
responding to low-order commensurabilities between them. They give origin to 
secondary resonances inside the primary 2/1 resonance, whose locations on the 
representative plane are indicated in Figure \E\ 

For some initial conditions, the amplitude (not the frequency) of one mode 
tends to zero, and the solution changes from quasi-periodic to periodic. Since 
we have two independent frequencies, we can also have two independent families 
of periodic motion inside the 2/1 resonance, which we refer to as a-family (the 
amplitude of the resonant mode is equal to zero) and Aw -family (the amplitude 
of the secular mode is equal to zero). On the dynamical map shown in Figure 
[5j the periodic families appear as continuous narrow white strips inside the gray 
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Figure 5. Dynamical map of the domain around the stable symmetric 
(0,0)-ACR (red dot) with e\ = 0.1 and m 2 /m 1 = 1.064. The initial 
value of G\ is fixed at and the initial value of Aw is fixed at (positive 
values on the e2-axis) or 180° (negative values on the e2-axis). 



tone domains of quasi-periodic motion. By definition, the intersection of two 
families gives us the position of the central ACR solution. 

Apart from the characteristic frequencies, the resonant and secular modes 
also differ in dynamical behavior. The resonant mode is always librating and the 
passage to circulation must occur at a separatrix, which gives origin to domains 
of highly chaotic motion, with possible disruption of the system. At variance, 
the secular mode is generally oscillating/circulating; only in the domains of very 
high eccentricities it is possible to find initial conditions which lead to a true 
secular resonance inside the 2/1 MMR (see Section [7.2. p . 

The domains of chaotic motion (and large instabilities) are always present on 
dynamical maps of the 2/1 resonance. The chaotic behavior is associated with the 
existence of separatrices between: (i) the 2/1 resonance region and the regions of 
near-resonant and purely secular motion; (ii) the regions of qualitatively distinct 
regimes of motion inside the 2/1 mean-motion resonance, and (Hi) secondary 
resonances inside the primary 2/1 resonance. The domains of highly chaotic 
motion are always shown in dark tones on the dynamical maps, while the hatched 
regions are regions of large-scale instabilities followed by disruption of the system 
within the time interval of each simulation (130,000 years). Finally, the domains 
of forbidden motion are filled by dark gray color in Figure [5j 
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Figure 6. Left panel: Dynamical map on the representative plane 
obtained for o\ fixed at and Aw at (positive values of the e2-axis) 
or 180° (negative values of the e2-axis). The position of the unsta- 
ble ACR is shown by a blue circle. Right panel: Dynamical map on 
the asymmetric representative plane obtained with initial values of the 
angular variables equal to those an asymmetric ACR: o\ = 45.4° and 
o"2 = —68.4°, for m^jmi = 0.3. The position of the stable ACR is 
shown by a red circle. 



7. Dynamics around the ACR solutions: Exterior resonance 

The dynamics associated to all asymmetric ACR is characteristic of exterior res- 
onances. Asymmetric stationary solutions arise for mass ratios mil 'mi lesser 
than 1.015 (see Figure [21 right panel). Moreover, for mi I mi < 0.36 and at high 
eccentricities, each asymmetric family bifurcates into two disconnected branches 
on the (ei, e2)-plane, referred to as branches A and B in our work. The rami- 
fication of the asymmetric families indicates the advent of the secular resonance 
inside the mean-motion resonance. As a consequence, the phase space of the 
exterior resonance shows a very complex picture marked by the presence of sev- 
eral distinct regimes of resonant and non-resonant motion, crossed by families of 
periodic orbits and separated by chaotic zones. 

7.1. Dynamical map around the one-branch asymmetric ACR 

One-branch asymmetric solutions are characteristic for the mass ratio values in 
the range from 0.36 to 1.015. The dynamics around a single asymmetric ACR is 
illustrated in Figure [6] on the two representative planes: the left graph is a plane 
analogous to one introduced for the interior resonance and shown in Figure [5j 
the right graph is a plane defined by the initial values of the resonant angles 
fixed exactly at those of the stable asymmetric stationary solutions. 

The plane defined by symmetric resonant angles does not allow us to fully 
describe the dynamics of the exterior resonance. The basic properties of the 
resonant phase space, such as the existence of the intersecting families of periodic 
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orbits and, thus, of the stable centers, are not seen on the resulting dynamical 
map. However, the analysis of this map allows us to detect the co-existence of 
the two regimes of motions named as interior resonance and exterior resonance 
for the same set of the constants m-zlmx-, AM and IC. Note the Aw-families of 
symmetric periodic solutions (with Aw = or 180°), which appear as narrow 
light strips inside the chaotic regions on the map in Figure [UJ left panel. 

To better represent the dynamics in the neighborhood of the asymmetric 
centers, we re-calculated the dynamical map for the same sets of constants of 
motion, but with the initial values of the resonant angles fixed exactly at those of 
the stable asymmetric stationary solutions (right graph in Figure [6]). Now we can 
observe the existence of the stable center at the intersection of the two periodic 
families (white strips), <r and Aw. The stable domain surrounding the center (in 
light tones), is composed of quasi-periodic motions with two independent modes: 
resonant libration around the cr-family and secular oscillation/circulation around 
the Aro-family. Generally, the global dynamics of the single exterior resonance 
is similar to the structure of the interior resonance presented in the previous 
section. However, there are two main differences: (i) the librating angle is o"2 
and the oscillating/circulating angle is o\ (and Aw), and (ii) the angles oscillate 
around values that are not simply or 180°. 




Figure 7. Smoothed planetary paths around one-branch asymmetric 
ACR in the (e\, Aw) (left panel) and (e2, 02) (right panel) polar co- 
ordinates calculated with e\ = 0.09 and mj/rai = 0.3 and the same 
energy. The smoothing has been done using a low-pass filter. 

The planetary motions around a stable asymmetric ACR are shown in Figure 
[3 in the polar coordinates (e\,Aw) on the left graph and in the polar coordinates 
(e2, 02) on the right graph. Only the secular components are presented; the 
fast resonant oscillations were eliminated using a low-pass filter. All planetary 
paths enclose the asymmetric ACR whose coordinates are e\ = 0.13, e2 = 0.26, 
a 2 = 284° and Aw = 108°. The orbits close to the ACR are presented by black 
curves: they show oscillations of both the resonant angle 02 and the secular angle 
Aw. For initial conditions far away from the stable center, the amplitudes of 
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the oscillations increase and the angle Aw starts to circulate, as shown by blue 
curve in Figure left panel. In this case, the curve separating the oscillatory and 
circulatory solutions represents a solution passing through the origin e\ = 0; we 
say that the secular angle Aw (and o~\) is in the oscillatory /circulatory regime 
of motion. 

At variance, the red curve in Figure [7J right panel, separating the oscilla- 
tions of the resonant angle a 2 around the asymmetric ACR from the oscillations 
around zero, is formed by solutions asymptotic to the unstable symmetric ACR. 
Hence this curve is a true separatrix and the oscillations of o~2 around the stable 
asymmetric ACR are true librations. 

Outside the separatrix, the resulting orbits (as an example, the external 
black curve in Figure [7J are horseshoe-like orbits encompassing both asymmetric 
centers and appearing to oscillate around zero (analogous of orbits known from 
the dynamics of the asteroidal 1:1 mean- motion resonance with Jupiter). On the 
dynamical map on the right panel in Figure [6l the domain of the horseshoe orbits 
is separated from the central resonance region by the blue curve; this curve is the 
locus of the initial conditions corresponding to the true separatrix corresponding 
to the red color path in Figure [7J The detailed analysis of the dynamics in its 
close vicinity reveals that the transition across the blue curve is topologically 
discontinuous, but the associated large-scale instabilities followed by disruption 
of the system are not observed on the dynamical map. 



7.2. Dynamical maps around the two-branch asymmetric ACR 

Two-branch asymmetric solutions are characteristic for mass ratio values lesser 
than 0.36. In this case, there exist two distinct stationary configurations de- 
scribed by different sets of orbital elements, but leading to the same value of 
total angular momentum. To better illustrate the bifurcation phenomenon, we 
analyze and compare planetary paths in the vicinity of the two-branch ACR: one 
belongs to the branch A and other to the branch B. This is done in Figure [HJ 
where, in order to simplify the visualization, we concentrate solely on the secular 
behavior of the system. The paths of both planets were smoothed using low- 
pass filtering (resonant oscillations were eliminated) and were plotted in polar 
coordinates (ei, Aw). 

Two centers are now clearly observed in Figure [HJ instead of one center as 
in the case of the one-branch asymmetric ACR (see Figure [JJ. The black curves 
surround the center associated to the branch A of the periodic solutions, while 
the red curves surround the center on the branch B. Among these paths, we 
find an infinite-period separatrix formed by solutions asymptotic to an unstable 
ACR (saddle point). Outside the separatrix, all paths follow similar patterns 
and are structurally stable. The existence of the separatrix, not present in the 
case of a single ACR, implies that there are motions near one of the two-branch 
asymmetric ACR that actually correspond to a true libration around the ACR 
and no longer a simple circulation/oscillation. This asymmetric ACR belongs to 
the branch B. In its neighborhood therefore both Aw and o~i show true librations 
around the equilibrium values. In other words, the two-branch asymmetric ACR 
generates, inside the 2/1 mean-motion resonance, a true resonance of the secular 
angle Aw = o~\ — o~i- This novel behavior is absent in the cases of symmetric 
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Figure 8. Smoothed planetary paths around two-branch asymmetric 
ACR in the (ei, Azu) polar coordinates calculated with e\ = 0.09 and 
m2/mi = 0.3 and the same energy. The smoothing has been done using 
a low-pass filter. 



and one-branch asymmetric ACR (except for very-high, close to 1, eccentricities 
of the inner planet). 

Using the same constants of motion, but two distinct sets of critical angles, 
we constructed the dynamical maps on the two asymmetric representative planes, 
shown in Figure [9] (top panels). The positions of the stable asymmetric centres 
on the corresponding dynamical maps are marked by red dots. The dynamical 
map of the branch A solution (left-top panel) displays qualitative similarities 
with the one-branch ACR map shown in Figure [SJ right-top panel, mainly the 
existence of two families of periodic motion intersecting at the center and a large 
zone of the quasi-periodic motion inside the exterior resonance regime. The main 
difference is the existence of the separatrix whose location is indicated on the 
left-top graph in Figure [9j 

In order to better visualize a separatrix feature, we present in the same figure 
(bottom panels) the neighborhoods of the stable asymmetric ACR, plotting the 
energy levels of the resonant Hamiltonian. A bifurcation phenomenon is clearly 
observed in the behavior of the Hamiltonian around the branch A ACR on the 
left-bottom panel, where the energy level in blue color presents a saddle-type 
structure. The continuation of the same level in the domain of the branch B 
ACR appears on the right-bottom graph as a blue curve. The fixed point of the 
centre (red dot) on the left-bottom panel appears as the energy level in red color 
on the right-bottom graph; that is, the branch B stationary solution is located 
at the global maximum of the resonant Hamiltonian. 
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Figure 9. Top: Dynamical maps of the two possible asymmetric so- 
lutions: one with u\ = 32.2° and 02 = —62.4° belongs to the branch 
A (left panel), and other, with o\ = 120.7° and a 2 = —69.1°, to the 
branch B (right panel) . The positions of the stable asymmetric centers 
on the corresponding dynamical maps are marked by red dots. Bottom: 
The levels of the resonant Hamiltonian around the branch A ACR (left 
panel) and the branch B ACR (right panel). 



The region surrounding the branch B ACR (right-top panel) is a domain of 
true secular resonance, when the passage from oscillation to circulation of the 
angle Aw is chaotic. In this region, both the resonant angle 02 and the secu- 
lar angle Aw librate around the corresponding asymmetric values (in the case 
shown in Figure [9j they are 120.7° and 189.1°, respectively). The equilibrium 
configurations of the planetary orbits are similar to those shown in Figure [4j the 
closest approaches occur when the central star and the two planets form a tri- 
angle configuration. Our numerical simulations have shown that the planetary 
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motions in this region are very stable, even for very high eccentricities of the 
orbits. 

Finally, on both maps in Figure [9] there exist regions of the horseshoe-like 
orbits at initial conditions situated outside the cyan curves. The horseshoe orbits 
are not centered on any single ACR, but actually contain all of them including 
the single (or double) stable asymmetric ACR and the unstable symmetric ACR. 
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